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ARTICLE INFO ABSTRACT 








In the present study a comprehensive thermodynamic model and degradation based optimization framework for 
energy management of anode supported solid oxide fuel cell (SOFC) stacks are carried out. The optimization 
framework determines optimum operating conditions to maximize system productivity (energy generation over 
system lifetime) considering degradation mechanisms. The main degradation mechanisms in anode supported 
SOFCs are nickel coarsening and oxidation. In this study, the optimum operating conditions regarding these 
degradation mechanisms to achieve maximum productivity at different target lifetimes are derived. The results 
show that target lifetime has a significant impact on system productivity and optimum operating temperature 
and current density. Furthermore, SOFC optimum operating conditions as a function of target lifetime are de- 
rived. To show the effectiveness of the developed framework, model outputs are compared with two other 
operating strategies; a base case strategy that optimizes system operating conditions without considering de- 
gradation mechanisms and a strategy based on Department of Energy’s (DOE) 2016 fuel cell report. Results 
illustrated that degradation based optimization is more beneficial for improving the entire performance in long- 
term operation. For instance, system productivity is 7.4% higher in comparison with DOE strategy during 
40,000 h operating lifetime. It is expected that the proposed methodology will lead to more rapid commercia- 
lization of SOFC technology. 


Keywords: 

Degradation based optimization 
Solid oxide fuel cell 
Productivity 

Operating conditions 

Target lifetime 





1. Introduction prediction approach for proton exchange membrane fuel cell perfor- 


mance deterioration is proposed based on a multi-physical aging model 


The solid oxide fuel cell (SOFC) is entirely solid-state and highly 
efficient. It has many advantages as a power generation device. For 
instance, it produces no noise during operation since it has no moving 
parts. In addition, the SOFC has no strict requirements on fuel gas 
composition. Unlike some other types of fuel cells that need pure hy- 
drogen, the SOFC can use variety of fuels, most of which are hydro- 
carbon-based fuels such as methane and propane. The other advantage 
is that since SOFC operates at high temperatures, usually between 
500 °C and 1000 °C, expensive platinum catalyst is not required. Also, 
the heat generated during operation can be used as a source of heat 
energy such as heating buildings [1]. However, due to issues caused by 
contaminants and degradation mechanisms over the long term, SOFCs 
are still not broadly implemented [2]. Thus, a fundamental under- 
standing of degradation mechanisms for long term analysis of SOFC 
system is needed. Several studies have been conducted in recent years 
that model degradation mechanisms and the effect of degradation on 
performance deterioration [3-5]. For instance, in [6], a novel 
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with particle filter approach. Optimizing operating conditions with 
considering performance deterioration leads to more accurate and re- 
liable results [7,8]. 

In this study, a degradation based optimization (DBO) framework is 
proposed [9], and optimum operating conditions are derived in order to 
advance the commercialization of SOFCs. 

Degradation based optimization is a model that accounts for system 
degradation mechanisms in the optimization procedure. The goal of the 
DBO is to derive the optimum operating conditions or design para- 
meters of the system to maximize or minimize a specific objective 
function. In energy engineering systems, the optimum operating con- 
ditions mostly result in minimizing the system total cost or maximizing 
the system energy production through system operating lifetime 
[10,11]. 

There is a significant body of research that has focused on SOFC 
optimization models [12]. In most studies the optimization strategies 
are different. In addition, they are optimizing either operating 
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Nomenclature 

A area, m° 

Co, gas-phase concentration of O2, mol m`?’ 
D mass diffusivity, m° s71 

E activation energy, kJ mol’ 

F Faraday’s constant, C equiv’ + 

H enthalpy, J mol? 


2 


i current density, Am 

I working current, A 

Kox rate constant for oxidation reaction, m? kg™ + s71 
K can nickel particles growth rate 

L length, m 

M molecular weight, kg mol” * 

n equivalent electron per mole of reactant, equiv mol * 
ni number of nickel particles, mol 

p pressure, Pa 

P dimensionless pressure 

pr site occupation probability 

R universal gas constant, J mol + K7 t 

r nickel particles radius 

Tox rate of oxidation, molkg~* s~! 

S° standard entropy, J mol ~t K7! 

T temperature, K 

TPB triple phase boundary, m m~? 

u decision variable 

v voltage, V 

Wni nickel particles volume 

x state variable 

X molar fraction of component in the mixture 
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b Ohmic resistance, Q 

n overpotential, V 

o electrical conductivity, Q7! m7! 

o porosity 

p tortuosity 

v diffusion volume of simple molecules, cm? 
a fraction of the reaction heat that is generated at the anode 
ô pre-exponential factor, Q7! m~? 

Oo bulk oxygen coverage, dimensionless 

y exponential activity parameter, kJ mol ~ + 
T perimeter 


Subscript and super script 


0 initial condition 
1 bipolar plate (interconnect) 
2 fuel channel 

3 anode 

4 electrolyte 

5 cathode 

6 air channel 

Act activation 

ch channel 

Conc concentration 
Ohm Ohmic 

Opt optimal 

X x direction 

y y direction 

Z z direction 








conditions or design parameters [13]. For instance, Feng et al. [14] 
developed a design optimization model to maximize system output 
power. In this study, the structure of a single tubular SOFC is considered 
as a decision variable. Results show that optimum structure leads 
system to operate at 18.2% higher output power in comparison with 
not-optimized structure design. Duhn et al. [15] proposed a novel op- 
timization model to determine the optimal geometric flow design. The 
modeling of the system is based on the computational fluid dynamic 
modeling and decision variables are the channels width and the area in 
front of the parallel channels. The optimal design maximizes flow 
uniformity index which is directly proportional to the maximum ob- 
tainable overall conversion. Hasanabadi et al. [16] investigated a model 
to optimize the microstructure design of an SOFC. In their study two- 
point correlation functions are used to realize the three dimensional 
porous microstructure of the SOFC. The results show that the optimi- 
zation procedure can be used as a robust tool to design the optimal 
microstructure. 

In [17] a multi-objective optimization model is investigated that 
optimizes operating conditions of a tri-generation system driven by 
SOFC. In this study the objectives of the optimization are to minimize 
system total product unit cost and maximize exergy efficiency, si- 
multaneously. The optimization algorithm is based on genetic algo- 
rithm and results indicate that optimal conditions are on the Pareto 
front. In [18], the operating conditions of a tri-generation system is 
optimized to make the system economically affordable while meeting 
the whole cooling demand and 50% of the total electricity demand. The 
model is optimized for a tri-generation system includes 50 kW tubular 
SOFC combined with heat recovery steam generator, combustion 
chamber and a chiller. Results show that optimal operating conditions 
lead system to be economically affordable and have an approximately 
net annual profit of $874,200. 

The design parameters and operating conditions can be optimized 
more accurately by considering degradation mechanisms in the 
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optimization procedure [19]. Xu et al. [20] developed an optimization 
model considering system degradation to optimize thermal stress per- 
formance of the SOFC. The aim of the model in their study is predicting 
the thermal stress of an anode-supported SOFC for different designs. A 
three dimensional thermal stress model is considered as the degradation 
mechanism. The results show that as the interconnect area increases, 
the thermal stress decreases. In addition, the thermal stress of the 
counter-flow and co-flow are similar. 

As is clear from reviewed literature, most SOFC degradation models 
are based on experimental data. However, in principle-based models, 
the degradation mechanisms rate dependency to operating conditions 
and SOFC material properties can be determined more accurately. 

As can be seen, in the case of SOFC optimization models, significant 
research has been done in the last decade [21-24]. However, there is 
not many studies that optimize system considering degradation 
[25-28]. The main contributions of the present study are as follow: 


e One of the main novelties of the proposed framework is the con- 
sideration of degradation model in the optimization procedure. The 
degradation mechanism can result in performance system long term 
productivity deterioration. The rate of degradation mechanism de- 
pends on system operating conditions. Effective control of de- 
gradation mechanism rate is a key element of the efficient operation 
of energy systems. The proposed model manages degradation me- 
chanism by optimizing system’s operating conditions in order to 
maximize lifetime productivity. 

e An improved degradation based optimization approach is used in 
order to address degradation mechanism effect on system pro- 
ductivity deterioration. This new approach results in more accuracy 
in long term simulation and more reliable and realistic optimal 
operating conditions. 

e The proposed framework produces optimal operating conditions of 
SOFC stacks to achieve maximum productivity over system target 
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lifetime. 

e The optimal operating conditions as a function of target lifetimes are 
derived in order to maximize system lifetime productivity. 

e In addition, a 150 kW SOFC power plant is modeled as a case study, 
and the optimal operating condition is derived. Results show the 
effectiveness of the proposed model and solution algorithm. 


For this purpose, a brief overview of the developed degradation 
based optimization framework is presented in Section 2. The formula- 
tion of SOFC degradation based optimization is presented in Section 3, 
and the sub modules are introduced. Results are presented in Section 4 
and finally the study is concluded with an overall discussion on the 
advantage of DBO framework development in Section 5. 


2. Generic problem statement 


Degradation mechanisms deteriorate system productivity over long 
term operation. On the other hand, degradation mechanisms rate is a 
function of operating conditions. Therefore, productivity deterioration 
over long term operation can be controlled by optimizing system op- 
erating conditions. As mentioned before, DBO model optimizes system 
operating conditions with considering degradation effects in the mod- 
eling. This model has two main components. Firstly, it studies the de- 
pendency of degradation mechanisms rate to the operating conditions. 
And secondly, it considers degradation mechanisms effect on system 
long term productivity (Fig. 1). 


Operating conditions 


Ko 


Degradation mechanisms 


| (2) 


System long term productivity 


Fig. 1. Two main parts of DBO model. 


In order to consider these two main components in the optimization 
model, system process, degradation and optimization modules should 
be developed. These three modules are integrated and solved con- 
currently. 

The process module is an Input/output (I/O) model of the system 
which formulates the dependency of system outputs such as perfor- 
mance to the system inputs like ambient and operating conditions. 

The degradation module can be developed based on the degradation 
mechanisms governing equations and available data. 

The optimization module, consists of an objective function and 
several constraints. The objective function shows the operator strategy 
of operation. In this study, maximizing power generation is the objec- 
tive of the operation optimization, and the constraints consider the 
operation and economical limitations of the system. 

The geometric feature of the studied SOFC is illustrated in Fig. 2, 
which includes external dimensions (Lx, Ly, and Lz) and individual 
layers’ thicknesses (Lj). 

The DBO model of a SOFC stack consists of process, degradation, 
and optimization modules. In the process module, system output power 
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Fig. 2. 3D configuration of a planar SOFC single cell. 
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Fig. 3. Flow diagram of a SOFC process module. 


as a function of operating conditions (cell temperature and current 
density) and design parameters is derived. The data flow diagram of the 
process model is presented in Fig. 3. 

The second module is degradation module that models the main 
degradation mechanisms including nickel coarsening and oxidation of 
SOFC stack. In this module, degradation rate of catalyst layer is cal- 
culated at different operating conditions. In addition, the effect of de- 
gradation mechanisms on stack design parameters deterioration such as 
triple phase boundary and conductivity is determined. 

System performance deterioration as a function of operating con- 
ditions can be calculated by integrating process and degradation 
modules. The data flow diagram of the integrated model is presented in 
Fig. 4. At each time step, degraded design parameters are derived from 
degradation module. The degraded design parameters are used as up- 
dated inputs of the modules for the next time step. 

The third module is the SOFC optimization module. In this module, 
operating conditions are optimized in order to maximize system life- 
time productivity. The objective function of optimization model is 
maximizing lifetime productivity and the decision variables are the 
operating conditions. Optimization model constraints are derived from 
process and degradation modules. The data flow diagram of SOFC DBO 
model (integration of process, degradation and optimization modules) 
is presented in Fig. 5. 
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Fig. 4. Flow diagram of SOFC integrated process and 
degradation modules. 


Next time step 


e Degraded catalyst layer 

e Degraded design parameters 
o Triple phase boundary 
o Electrical conductivity 





Fig. 5. Flow diagram of SOFC DBO model. 
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3. Degradation based optimization model 
3.1. SOFC process module 


The flow diagram of process module is presented in Fig. 6. Inputs 
and operating conditions are applied to the voltage calculation sub- 
module and system output power is derived based on the governing 
equations presented in this section. 

The main governing equation is the actual operating potential 
evaluation (V), (Eq. (1)). This term is derived from the open circuit 
potential (Vocy) by considering occurred losses include: Ohmic Mam)» 
concentration (onc) and activation (7,.,) over potentials [24,25]. 


V (t) = Vocv—( 2 Nohmj (t) + Neonc,3 (t) + Neonc,5 + Nact,3 (t) + Nact,s) 
J 


j = 1,3,4,5 (1) 


The open circuit potential is estimated by Eq. (2) [24,25]. 
Vocv 


P 
—(Hio-Hin—>Ho,) + T x (Sf,o-Sfn—75°O2)-RTIn | tae 
H (t 02 


2F 
(2) 
The Ohmic overpotential is calculated based on the Ohmic law. The 


Ohmic overpotential of interconnect and electrolyte are calculated as 
Eq. (3) [24,25]. 
Nonm j = Ifj = 14 (3) 


This loss depends on the electrical resistances of interconnect (£) 
and electrolyte (6,) which are calculated from Eqs. (4) and (5) [24,25]. 








Lı Lz Le 
p= taoa 
LyLz0 (1 + Ach) Lt Lyo (1 + Nen)LiLz0 (4) 
L 
By = 
LyLz04 (5) 


Moreover, the ohmic overpotential of cathode and anode are cal- 
culated as Eq. (6) [24,25]. 


I(Aj,opt (t) — LyLz)Lj,opt (t) I(Lj — Lj,opt (t)) 


Moimi = FA opi()iylz = ga t Tyla- LO 
j= 3,5 (6) 
r=------------ ; 
, Operating conditions ! 
| «Current density i 
i e Temperature 


Process module 


Open circuit potential 
Eq. (2) 
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In the Eq. (6), the area of the active TPB regions (Aj opt) and average 
diameter of electrodes parallel ports (Dp;) are calculated as follow 
[24,25]. 


Aj opt (£) = (1 F 2L; opt (t)/Dpj)Lyz (7) 


(8) 


The concentration over potentials in the cathode and anode layers 
are given by Eqs. (9) and (10), respectively [24,25]. 


RT Po, 
as ee ae 
4F \ Po, rps 


| Pr, Pro, rep (t) 


Prpo,rep(t) Pro 


Dp; = Dwj®;/(A—-G;) 


Nconc,5 = 


(9) 


RT, 


Neonc,3 (t) = 4F 


(10) 

The required partial pressure of water (Py,0,rpg), partial pressure of 
hydrogen (Px, rpg), partial pressure of oxygen (Po,,rpg), mass diffusivity 
at anode (D3) and cathode side (D;) for calculating concentration 
overpotentials are obtained from Eqs. (11-15) [24,25]. 


RTP, L;I 





Pryo,rep(t) = PpXmo + 





FDD, LL B 
RTh, L 
Prp,rep(t) = Pe Xma 
2 p 2FD3(t)¢, LyLz a 
A i x ) RTġ,.L;I 
=p. —(ņ. — ip )EXP| -n n yr 
O2,TPB = Pair —Pair O2 Pair P 4FD5 5 Pai Ly Lz (13) 
1.43 x 10-7T"5 (My, + Mm o)2 
Pr (Mn, Mmo): (vå E vko) (14) 
1.43 x 107-7775 (Mo, + My,)2 
Pair (MoM) (06 + v8] (15) 


The last over potential is activation which is given by Eq. (16), 
where ig; stand for the exchange current densities and are calculated 
from Eq. (17) [24,25]. 


RT [ 
BE ante i 
anF i 





Nactj (t) = 
(16) 





l 
i 
l 
e Output power - 
Eq. (18) i 

1 


Fig. 6. Detailed data flow diagram of a SOFC process module. 
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RT E; 
io; (t) = —kexp| —— j = 3,5 
of) = pF r( z) l (17) 
Based on the calculated actual operating potential (Eq. (1)), the 
electrical power density of SOFC is obtained by the multiplication of the 
actual operating potential and current density as represented in Eq. 


(18). 
Power(t) = V(t) xi (18) 


3.2. SOFC degradation module 


As system operates, different degradation mechanisms occur. The 
case being studied is an anode supported SOFC stack with stationary 
application. The main degradation mechanisms in anode supported 
SOFCs are nickel coarsening and oxidation which occur at the anode 
side. Nickel phase in the yttria-stabilized zirconia (YSZ) structure will 
be reformed in unsuitable operating conditions. Nickel particle size 
increases in this undesired reformation mechanism. Moreover, nickel 
will be oxidized in the presence of oxygen, water and oxygen ions at 
anode side. The increase in nickel particle size and nickel oxidation, 
diminish the contact of reactants within the Ni phase. Therefore, the 
triple phase boundary and the electronic conductivity are degraded 
[27-30]. 

Based on the degradation reactions, the flow diagram of degrada- 
tion module is presented in Fig. 7. The model is run at each time step 
and parameters are updated for the subsequent time step. At the be- 
ginning of each time step, inputs and operating conditions are applied 
to the governing equations of degradation module and nickel coar- 
sening and oxidation rate are derived. The rate of these two degrada- 
tion mechanisms determine remaining catalyst active surface area. As a 
result, degraded design parameters including triple phase boundary and 
electrical conductivity are calculated. These values are the updated 


' Operating conditions 


l i 
l l | 
I | 
i | e Current density T 
e Temperature | 


Degradation module 


Nickel coarsening 
Eqs. (19, 20) 


Nickel oxidation 
Eqs: (21,23) 


Degradation module output 


Degraded catalyst layer 
Eq. (24) 


Triple phase boundary 
Eqs. (25-27) 


Egs. (25, 26, 28) 
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design parameters, and are the inputs of the next time step. 

As can be seen in Fig. 7, nickel coarsening and nickel oxidation are 
the two main sub-modules of degradation module. Tanasini et al. [31] 
developed a nickel coarsening model, in which particle coarsening 
obeys charging capacitor law as Eq. (19). 


R(t) = (R(max)—R (0))[1—-exp(—Ke,capt) | + R(0) (19) 


R(t) is nickel radius at time t. R (max) and R (0) are maximum and initial 
radius of nickel particles. And, ks,cap is the capacitor constant. According 
to mass balance of nickel particles on catalyst layer, number of re- 
maining nickel particles can be calculated as Eq. (20). 


3 
i= (30) x ni(0) 


R(t) (20) 


Moreover, nickel can be oxidized by the presence of oxygen or water 
vapor at anode side [26]. 


1 
Ni + —O, © NiO 
2°? (c) 


Ni + H, Oe NiO + H, (d) 


To evaluate the oxidation rate of each reaction, mole fraction of O, and 
H,0 at anode side should be determined. The amount of H,O depends 
on input fuel and oxidant. However, oxygen can exist at the anode side 
because of leakage or water electrolysis [27,28]. The nickel oxidation 
rate is calculated as Eq. (21). 


dni(t) 
ox = = Kox 1-9, i 
fox = aa (21) 


Where k,, is the reaction rate coefficient and calculated as Eq. (22). 


Box Y 
Kox = k’ = =p 
exp( = exp RT 








(22) 


Fig. 7. Detailed data flow diagram of a SOFC de- 
gradation module. 


Next time step 
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The remaining nickel particles is calculated as Eq. (23). 


ni (t) = ni,(t)— J. Fox» dt (23) 


In this equation, ni,(t) is the remained nickel particles with con- 
sidering coarsening mechanism. The available active surface area of 
catalyst layer is corresponding to available nickel particle surface area 


that is calculated as Eq. (24) 
A(t) = ka X ni (t) x 47R(t)? (24) 


In addition, nickel coarsening and oxidation decrease occupation 
probability and as a result affect the available nickel particles. 


pr(t) _ A(t) 
Pr, A(0) oa) 
Tri (t) = o| PO) + 3.6ni(t)° 
pr (t) (26) 


And finally, triple phase boundary and electrical conductivity of 
SOFC are calculated as Eqs. (27), (28), respectively [29]. 


TPB(t) = 2[(a—28)r (t) ]ti(T) (27) 


Oop (t) = oo (tan@) [pr (t)—pr.]|'° (28) 


In these equations, r(t), ty;(t) and pr(t) are degraded over the op- 
erating lifetime and as a result triple phase boundary and electrical 
conductivity deteriorate over time. These deteriorated parameters are 
the input of process and optimization modules. 


3.3. SOFC optimization module 


The optimization objective is to maximize total power generation 
through system lifetime (productivity) with considering degradation 
mechanisms. 

The fuel cell stack is the most expensive component with a higher 
failure rate in comparison with other components of a SOFC power 
plant [32]. Therefore, a complete replacement of the stack is considered 
as system target lifetime and it is assumed that other components can be 
maintained or replaced with a low cost during stack operation [32,33]. 
In addition, it is assumed that the operating condition of the SOFC 
power plant components are independent. 

As a result, in this study, the main part of the power plant which is 
SOFC stack is studied and optimal operating conditions of the stack is 
derived. 

Decision variables of the optimization model are the operating 
conditions of the stack which are cell temperature and current density. 
The optimization module flow diagram is presented in Fig. 8. 

The objective function is presented as Eq. (29). 


t+T 
max J [u(t),x(t)] = J Power (t)dt (29) 


Power (t) is the SOFC output power which here is as Eq. (30). 


Power (t) = V(t) xi xX A(t) (30) 


V(t) and A(t) are the actual operating potential and active surface area 
which are derived based on process and degradation modules. The 
optimization constraints set includes restrictions on the system opera- 
tion. For instance, the operating temperature and current density 
ranges are considered. Moreover, it is assumed that SOFC operates up to 
half of the nickel degraded [34]. 

To solve the proposed degradation based optimization model, the 
genetic algorithm (GA) method is used [35]. The GA algorithm starts 
from a population of randomly generated individuals (operating con- 
ditions) and occurs in generations. In each generation, the optimization 
output (productivity) is calculated, and the fitness of every individual in 
the population is evaluated. Multiple individuals are selected from the 
current population (based on their fitness), and modified to form a new 
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population. The new population is used as an initial value in the next 
iteration of the algorithm. The loop is repeated until the required 
number of iterations yield the same optimal results. 


4. Results and discussion 
4.1. Case study 


As a case study, a 150 kW SOFC is selected. Properties and constant 
values used in the simulations are referenced from Wen et al. [25] re- 
search and are presented in Table 1. 


4.2. Validation 


The model should be validated before it is applied to the case study. 
Two steps of validation are performed to illustrate that results are re- 
liable. In the first step, to validate the SOFC process module, compar- 
ison between model outputs with results published by Wen et al. [25] is 
performed. Fig. 9 illustrates that under the same SOFC geometric di- 
mensions and operating condition (1073 K), the general shape of two 
polarization curves are similar. 

According to the results, the absolute error varies between O and 
5%. If the purpose of the simulation model was to have an as accurate 
model as possible, this might have been problematic, but in this case 
our focus is to use the model for operation optimization, and deriving 
an accurate parameters trend is sufficient, which is accurately achieved. 

In the second step, the integration of SOFC degradation and process 
modules is validated. The voltage degradation rate at different current 
densities is used to validate this integrated model. As mentioned before, 
the main degradation mechanisms in anode supported SOFCs occur at 
anode side. Therefore, the integrated model results are compared with 
Muller et al. [36] experimental research that studies anode degradation 
mechanisms. In [36], the efficiency and long term stability of Ni-YSZ 
anodes as a function of current density are studied for various SOFC 
single cells. Tables 2 and 3, compare experimental data gathered by 
Muller et al. [36] with modeling results at two different operating 
current densities representing the low and high ranges of operations. 

Comparison at low current density range, shows that the voltage 
degradation rate is slightly higher in the simulation than what is 
measured. 

In addition, the absolute error at higher current densities is greater 
due to the existence of other degradation mechanisms such as cathode 
degradations which are ignored in the modeling. 

Improvements in the agreement between the simulations and ex- 
perimental data can be obtained by modeling all degradation me- 
chanisms such as cathode degradation mechanisms. However, the main 
objective of the present study is to capture degradation profile over 
time which is achieved with enough accuracy. 


4.3. Optimal operating points 


Solid oxide fuel cell operating conditions (temperature and current 
density) have a great impacts on power deterioration over time. SOFC 
can operate at operating temperature range of 873-1123°K [37]. And, 
the current density ranges of 1-10 A/cm? [38]. 

In order to maximize system power generation through lifetime, the 
optimum operating conditions should be determined. In addition, 
system target lifetime affect the optimum values. In this study, optimum 
operating conditions are determined based on the developed DBO 
model for different target lifetimes. 

The optimized operating conditions and power generation through 
time according to the three target lifetimes are shown in Fig. 10. When 
the target operating time is 10,000 h, the initial power is 136.5 kW at 
the optimum temperature and current density of 1065.2 K and 2.69 A/ 
cmĉ?, respectively. As target lifetime increases, the initial power and 
optimal operating conditions decrease. For instance, at 40,000 h target 
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Degradation module output 
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Max Productivity 
Eqs. (29) 





Constraints 





Output power= f ( Design parameters , Decision variables) 


Target lifetime= Const. 


Degradation rate = f ( Design parameters , Decision variables) 
Lower bound < Decision variables < Upper bound 


Eqs. (1-28, 30) 
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Fig. 8. Detailed data flow diagram of SOFC opti- 
mization module. 
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Table 2 
Comparison between experimental data and modeling results at low current density. 
Current density 


Anode voltage degradation rate (V/ Differences (V/ 


(A/cm?) 1000 h) 1000 h) 
Experimental data Modeling results 
0.2 0.032 0.039 + 0.007 
Table 3 


Comparison between experimental data and modeling results at high current density. 


Current density Anode voltage degradation rate (V/ Differences (V/ 








(A/cm?) 1000 h) 1000 h) 
Experimental data Modeling results 
0.52 0.057 0.047 — 0.01 
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Fig. 11. Optimized initial power at different target lifetimes. 


operating time, the initial power is 114.7 kW at optimum temperature 
and current density of 1039.6 K and 2.11 A/cm”, respectively. These 
results suggest that the operating conditions should be decreased for 
long term operation but increased if the expected lifetime is shorter. In 
other words, increasing temperature and current density are more 
beneficial for improving the entire performance in short-term opera- 
tion. 

Fig. 11 presents the optimum initial power at different target life- 
times. As is clear, at higher target lifetime, system operates at lower 
optimum initial power because of lower degradation rate in this range 
of power. 


89 


Energy Conversion and Management 158 (2018) 81-91 


1070 2.8 


Optimized Current density = Se-10t? - 4e-5t + 3.05 





N P 
1065 g 
Z Q 
= 242 
É 1060 > 
z z 
Ə 1055 3 
E 2 § 
1050 p 
: 3 
= z 
a 1045 =| 
Š 16 
2. 
1040 Q 
Optimized temperature = 2e-8t? -1.8e-3t + 1080.5 
1035 1.2 
0 10,000 20,000 30,000 40,000 50,000 


Target lifetime (hr) 
+.. Temperature — Current density 


Fig. 12. Optimum operating conditions as a function of target lifetime. 
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Fig. 13. Catalyst layer degradation and power deterioration rates as a function of target 
lifetime. 


4.4. Optimal operating functions 


The DBO model is solved for different target lifetimes and optimum 
temperature and current density are derived for each target lifetime. 
Then, discrete optimal points are fit to two curves by regression analysis 
as presented in Fig. 12. These curves present optimal operating condi- 
tions (temperature and current density) as a function of target lifetime. 

As presented in Fig. 12, by increasing target lifetime, the decrease in 
the optimal current density is slower, whereas the optimal temperature 
decreases consistently. Based on the derived curves, the dependence 
function of temperature and current density to the target lifetime are 
determined and presented in Fig. 12. 

Fig. 13 shows the optimal output power deterioration and the de- 
gradation rate of catalyst layer active surface area as a function of 
target lifetime. When the target operation life time is 10,000h, the 
power deterioration rate per 1000h is around 0.9% when operating at 
optimum temperature and current density of 1065.2 K and 2.69 A/cm? 
(derived from Fig. 12). Fig. 12 illustrates that by increasing target 
lifetime, optimal operating temperature and current density decrease 
rapidly. Based on the Eqs. (21), (22), nickel oxidation rates is lower at 
low temperature and current density ranges. The low nickel degrada- 
tion rate results in lower degradation rate of catalyst layer (Eqs. (23), 
(24)), and power deterioration consequently. This phenomenon can be 
seen in Fig. 13. As target lifetime increases, catalyst layer degradation 
rate and consequently power deterioration falls. When the target op- 
erating time is 25,000 h, the performance deterioration is around 0.32% 
per 1000h at optimum temperature and current density of 1048.2 K 
and 2.29 A/cm”. These results suggest that the operating temperature 
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Fig. 14. Comparison of optimal productivity and operating conditions at different target 
lifetimes for three operation strategies. 


and current density should be decreased to achieve lower performance 
deterioration for long operation but increased if the required lifetime is 
shorter. This trend is the same for catalyst layer degradation rate as 
well. As shown in Fig. 13, in accordance with increasing target lifetime, 
catalyst layer degradation rate decreases and vice versa. It means that 
increasing the operating temperature and current density is more 
beneficial for improving the entire performance for brief operation 
because nickel degradation rate is high in this range of operating con- 
ditions and system will not work efficiently for a long duration. 


4.5. Comparison 


Fig. 14 illustrated the comparison between three operation strate- 
gies for the SOFC. Base case strategy optimizes stack operating condi- 
tions without considering degradation mechanisms in the optimization 
procedure. DOE strategy is based on Fuel Cell Multi-Year Research, 
Development, and Demonstration annual report [39]. According to this 
report, the performance degradation rate of SOFC stack with stationary 
application is 2% in 2016, [39]. In addition, the operating lifetime of 
the cell is between 12,000-70,000 h depending on operating conditions 
of the system, [39]. 

The productivity of the SOFC stack operating at different strategies 
is calculated for different target lifetimes of 10,000, 25,000, and 
40,000 h. As can be seen from the figure, DBO model optimum values 
lead to a more efficient and durable operation of the system than the 
other operation strategies. Based on the results, at 10,000h target 
lifetime, productivities of DBO and DOE strategies are nearly the same 
and working in DBO operating strategy leads to 27.8% higher power 
generation in comparison with base case strategy. 

As target lifetime increases, the difference between productivities 
increases. The base case strategy productivity is much lower than two 
other strategies due to ignoring degradation mechanisms in optimiza- 
tion model. It should be noted that the criterion for cell death is the 
degradation of half of the nickel. Therefore, system fails after 11,900h 
at base case operating strategy and energy generation for both 25,000 
and 40,000 h operating time are the same. 

However, DOE strategy considers degradation mechanisms in 
system operation scheduling. This strategy doesn’t optimize operating 
conditions in long term operation and as a result has a lower pro- 
ductivity in comparison with DBO strategy. 


5. Conclusion 


The performance and degradation of the SOFCs are strongly de- 
pendent on the operating conditions (temperature and current density). 
Thus, it is very important to determine the proper operating conditions 
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to ensure efficient long term operation. The purpose of this paper is to 
predict the performance and durability of anode supported SOFC and to 
find the optimal operating conditions for the target operating time. To 
achieve these objectives, firstly theoretical process module is developed 
to estimate cell performance. Secondly, degradation module is devel- 
oped based on the main degradation mechanisms in the anode sup- 
ported SOFC and the effect of operating conditions on anode catalyst 
degradation is determined. Moreover, by considering process module in 
correlation with degradation module, the effect of anode catalyst de- 
gradation on performance deterioration is also derived. The results 
from all developed models corresponded to the experimental data very 
well. Finally, the optimum operating temperature and current density 
to satisfy the desired lifetime are found by applying optimization 
module to the developed process and degradation modules. 

Based on the developed DBO model, the optimum temperature and 
current density as a function of target lifetime are derived to achieve 
maximum productivity. As a result, at any given target lifetime, the 
plant operator is enable to implement the optimal operating conditions 
which maximize system productivity and profitability. Furthermore, to 
demonstrate the value of DBO framework, system productivity oper- 
ating at DBO strategy is compared with two other operating strategies. 
Results demonstrate that the DBO model outputs lead to a more effi- 
cient and durable operation of the system compared with the other 
operation strategies. The suggested methodology is also valuable for 
commercialization of the anode supported SOFCs. 
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